arXiv:1505.01768vl [physics.flu-dyn] 7 May 2015 


Geometric pumping in autophoretic channels 

Sebastien Michelin,^’[^ Thomas D. Montenegro-Johnson,^ 

Gabriele De Canio,^ Nicolas Lobato-Dauzier/ and Eric Lauga^ 

^ LadHyX - Departement de Mecanique, Ecole Polytechmque - CNRS, 91128 Palaiseau, France 
^Department of Applied Mathematics and Theoretical Physics, 

University of Cambridqe, Cambridqe CBS OWA, United Kinqdom 
(Dated: May 8, 2015) 

Many microfluidic devices use macroscopic pressure differentials to overcome viscous friction and 
generate flows in microchannels. In this work, we investigate how the chemical and geometric 
properties of the channel walls can drive a net flow by exploiting the autophoretic slip flows induced 
along active walls by local concentration gradients of a solute species. We show that chemical 
patterning of the wall is not required to generate and control a net flux within the channel, rather 
channel geometry alone is sufficient. Using numerical simulations, we determine how geometric 
characteristics of the wall influence channel flow rate, and confirm our results analytically in the 
asymptotic limit of lubrication theory. 


I. INTRODUCTION 

Controlled flow manipulation at the micro- or nano¬ 
scale is at the heart of recent developments in microflu¬ 
idics, including many applications in the field of biolog¬ 
ical analysis and screening [T]. Generating and control¬ 
ling a flow within the confined environment of a mi¬ 
crofluidic channel requires an external forcing to over¬ 
come the viscous stress on the walls. In synthetic mi¬ 
cro and nanofluidic systems, this is usually achieved ei¬ 
ther mechanically, by applying a pressure difference be¬ 
tween the inlet and outlet of the domain, or through elec¬ 
trokinetics/electroosmosis, where the flow results from an 
externally-imposed electric field within the channel [ng. 

However, many biological systems rely on stresses lo¬ 
calized at boundaries in order to drive flow, rather than 
on a global macroscopic forcing. For example, micro¬ 
scopic cilia on the lung epithelium induce a directional 
flow of mucus through their coordinated beating, acting 
as a pump [6]. Similar microscale flow forcing at the 
wall also plays an essential role in the early stages of em¬ 
bryo development [7] or in the reproduction of mammals, 
where cilia-driven flow is responsible for the migration of 
the ovum down the female reproductive tract [8]. 

In a dual process, cilia-driven flows play an essential 
role in the self-propulsion of micro-organisms such as 
Paramecium [5]; the flow generated by the beating of 
cilia anchored on the wall of a moving cell is responsible 
for its locomotion. For both swimming and pumping, the 
coordination of neighboring cilia into metachronal wave 
patterns has been proven essential to achieving maxi¬ 
mum flow rate/swimming speed with a minimum ener¬ 
getic cost pnuid] . 

Several attempts have been made to reproduce ciliary 
pumping in the lab through the fabrication of artificial 
actuated cilia [IlHIH]- All of them rely, however, on the 


* Electronic address: sebastien.michelin@ladhyx.polytechnique.fr 


application of a global electromagnetic forcing field, and 
generating efficient pumping would require the applica¬ 
tion of phase-shifted forcing on neighboring cilia [niiio]. 
This constraint, as well as the manufacturing process of 
microscopic cilia, poses important challenges to minia¬ 
turization. 

Phoretic mechanisms, namely the ability to generate 
fluid motion near a boundary under the effect of an ex¬ 
ternal field gradient, represent an alternative route for 
both pumping and swimming systems that require no 
moving parts. These mechanisms arise from the inter¬ 
action of rigid boundaries with neutral or ionic solute 
species in their immediate environment, and are known 
to generate the migration of passive particles in external 
gradients m- Phoretic motion has recently received re¬ 
newed attention in the context of artificial self-propelled 
systems. Such artificial swimmers generate the field gra¬ 
dients required for propulsion themselves, for example 
through chemical reactions catalyzed at their surface, 
and thus do not rely on any external forcing to achieve 
propulsion [22II25j . These systems combine two proper¬ 
ties: (i) an activity, i.e. the ability to release/consume 
solute species or thermal energy at their surface, and (ii) 
a phoretic mobility, i.e. the ability to create a slip velocity 
at the boundary from a local tangential gradient of so¬ 
lute concentration (diffusiophoresis), temperature (ther¬ 
mophoresis) or electric potential (electrophoresis). Re¬ 
cently, autophoretic systems have also been considered 
for generating micro-rotors that rotate without the ap¬ 
plication of external torques [27] . 

In order to generate the surface gradients necessary to 
their self-propulsion, autophoretic particles must break 
spatial symmetry. A similar requirement exists for au¬ 
tophoretic pumps. This symmetry-breaking may be 
achieved for self-propelled particles following three main 
routes: (i) chemical asymmetry, i.e. patterning the parti¬ 
cle surface with active and passive sites (e.g. Janus par¬ 
ticles) [23l [24l [28] , (ii) spontaneous symmetry-breaking 
resulting from the advection of the field responsible for 
the phoretic response by the flow itself diisn] and (iii) 
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geometric asymmetry mM- 

In this work, we use a combination of theoretical anal¬ 
ysis and numerical simulations to investigate whether, 
and how, autophoretic mechanisms and geometric asym¬ 
metry can generate a net flow in a microfluidic channel 
without imposing any external mechanical forcing, elec¬ 
tromagnetic forcing or chemical patterning. For simplic¬ 
ity we focus on the specific case of diffusiophoresis, where 
slip velocities are generated at the wall from tangential 
gradients in the concentration of a solute released from 
one of the channel walls into the fluid. Because of the 
similarity between the different phoretic mechanisms, it 
is expected that the results of the present contribution 
may easily be generalized to thermo- or electrophoretic 
systems. Specifically, we follow the classical continuum 
framework of self-diffusiophoresis [231 155H55] . and con¬ 
sider how a left-right asymmetry in the wall shape can 
generate a net flow within the channel which hence acts 
as a microfluidic pump. 

The paper is organized as follows. Section [^ summa¬ 
rizes this continuum framework for the case of the flow 
within an asymmetric channel, and presents the numer¬ 
ical methods used in this work. Section m shows how 
the wall geometric characteristics determine the net flow 
within the channel. The results are then confirmed ana¬ 
lytically in Sec. jlV] using lubrication theory in the long- 
wavelength limit, and conclusions and perspectives are 
presented in Sec. 

II. PROBLEM FORMULATION 
A. Diffusiophoretic channel 

We consider a two-dimensional channel of mean height 
H, bounded by a flat bottom boundary (y = 0) and a 
top wall with a periodic non-flat profile, y = h{x) (see 
illustration in Figure]^. In the channel gap, filled with a 
Newtonian fluid of dynamic viscosity fi and density p, a 
solute species of local concentration C'(x) with molecular 
diffusivity D is present and interacts with the channel 
walls through a short range potential. When the typi¬ 
cal thickness of this interaction region is much smaller 
than the other length scales of the problem (namely the 
channel gap and the wavelength), the interaction of the 
wall with a local solute gradient generates an effective 
slip velocity at the wall [211135] 

Usiip = MV||C, (I) 

where Vy = (1 — nn) • V is the tangential component of 
the gradient to the surface of local normal n and M, the 
phoretic mobility, is a property of the solute-wall inter¬ 
action which may be positive or negative depending on 
the repulsive or attractive nature of that interaction |2Ij . 
The chemical properties of the channel walls are also 
characterized by a chemical activity, i.e. the ability to 
create or consume the solute species. Here we consider 


a simple fixed-flux model, for which the activity of the 
wall is given by a fixed flux of solute per unit area A, 
counted positively (resp. negatively) when solute is re¬ 
leased (resp. absorbed) 

DnVC = -A. (2) 

In the case of self-diffusiophoretic propulsion, loco¬ 
motion is often achieved through inhomogeneity in the 
chemical treatment of the particle [221424] . Recent work 
has shown that geometric asymmetry of chemically- 
homogeneous particle alone is in fact sufficient to en¬ 
sure locomotion ED [32]. Here we investigate a similar 
question, namely the possibility of obtaining a net flow 
from chemically-homogeneous channel walls using shape 
asymmetry. We thus assume that the top corrugated wall 
has homogeneous mobility M and activity A. To ensure 
the existence of a steady state solution, the concentra¬ 
tion of the solute on the bottom wall is assumed to be 
fixed {C = Co). Consequently, the fluid velocity on that 
wall satisfies the no-slip boundary condition. By study¬ 
ing the relative concentration of solute to that on the 
bottom wall, we can assume without loss of generality 
that Co = 0. 

The phoretic slip velocity generated at the top wall 
by the wall-solute interaction drives a flow within the 
channel. When viscous effects dominate inertia (namely, 
when the Reynolds number Re = pUH/p, is small, with 
U = \AM\/D the typical phoretic velocity), the flow sat¬ 
isfies the incompressible Stokes equations 

= Vp, V • u = 0, (3) 

for the velocity and pressure fields, u and p respectively. 
Solute molecules diffuse within the channel, and in gen¬ 
eral can also be advected by the phoretic flows. How¬ 
ever, when diffusive effects dominate (i.e. when the Peclet 
number, Pe = UH/D, is small), the solute dynamics is 
completely decoupled from the flow, and the solute con¬ 
centration satisfies Laplace’s equation 

V^C = 0. (4) 

Equations pi"®, together with the boundary condi¬ 
tions Eqs. Q-® applied at the top wall and the inert 
boundary conditions C = 0 and u = 0 at the bottom 
wall, form a closed set of equations that can be solved 
successively for the solute concentration C and velocity 
field u= {u,v). From these results, the net flow rate 
within the channel, Q, can be computed as 

rh{x) 

Q= uix,y)dy, ( 5 ) 

Jo 

which is independent of x because the flow is incompress¬ 
ible. 

B. Asymmetric channel 

The shape of the channel is a periodic function of x 
characterized by a wavenumber k = 27r/L. A sinusoidal 
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FIG. 1: Asymmetric phoretic channel. The top wall is charac¬ 
terized by constant chemical activity A and mobility M. The 
bottom wall maintains a fixed concentration and thus flow 
satisfies no-slip there. The example shown here corresponds 
to 7 = 7r/4, a/H = 1/2 and L/H = 27r with the asymmetric 
shape of the wall given in Eq. in¬ 


wall will generate a perfectly left-right-symmetric con¬ 
centration distribution and flow pattern, leading to no 
net flow along the channel. In the following, we focus on 
a subset of asymmetric wall shapes, essentially smoothed 
ratchets, that are formally obtained by mathematically 
shearing the symmetric sinusoidal profile. The top wall 
is described in parametric form by 

a;(s) = s —sin As, y{s) = H + asinks, ( 6 ) 
k 

the non-dimensional asymmetry parameter, 7 , deter¬ 
mines the asymmetry of the profile, and H and a are 
the mean channel width and the amplitude of the width 
fluctuations, respectively (see Figure [^. 

Hereafter, the problem is non-dimensionalized using 
1/fc as characteristic length, U as characteristic veloc¬ 
ity, and \A\/Dk as characteristic concentration fluctua¬ 
tion. While A and M are both signed quantities, after 
nondimensionalisation we focus on the case A — M = 1; 
changing the sign of either A or M simply reverses 
the slip velocity forcing and flow rate without chang¬ 
ing its magnitude. The problem is now completely spec¬ 
ified by three geometrical quantities, namely the non- 
dimensional mean gap, H* = kH, the corrugation am¬ 
plitude, a* = ka, and the asymmetry parameter, 7 . 


C. Numerical method 


function for Laplace’s equation, Eq. Q, is given by 

1 °° 

^ \n[{x-^ + 2mr)'^+ {y-T])^] 

n— — oo 

= ^ ln[ 2 (cosh( 2 / - y) - cos(x - C))]- (7) 

Assuming that the channel walls are smooth, the concen¬ 
tration at a point {x,y) on one of the walls can then be 
computed using the boundary integral formulation [36] 


\c[x,y) = f 

^ Jon 






( 8 ) 


The upper and lower boundaries of the domain are 
discretized into 200 straight-line segments, and C and 
dC/dn are assumed constant over each element. For 
elements on the bottom (resp. top) boundary, (7 = 0 
(resp. dC/dn = — 1) is enforced at the midpoint of each 
segment. This reduces the boundary integral equation 
0 to a dense matrix system for the solution vector con¬ 
taining the unknown dC /dn on the lower boundary and 
C on the upper boundary. The free-space (i.e. singu¬ 
lar) component of the Green’s function is isolated and 
integrated analytically, and all non-singular element in¬ 
tegrals are computed with a 16-point Gaussian quadra¬ 
ture. In order to compute the fluid flow and flow rate in 
the channel, only the boundary concentration of solute 
(and not its bulk distribution) is needed, which makes the 
boundary element method particularly suitable for this 
problem. The numerical code was validated against ana¬ 
lytical solutions for diffusion in a channel with nontrivial 
boundary conditions and domain geometry, achieving a 
relative error of at worst 0.004%. 

For Stokes flow, the dimensionless boundary integral 
equation for boundary force density, f, is given by 

Wi(x) = ^ Ohio 

- Tyfc(x-|)uj(|)nfc(|)]ds(|). (9) 

For X = X — ^ and r = |x|, the two-dimensional, 27r- 
periodic Green’s functions for Stokes flow are 


S = Y, I In - 


T = 


n— — oo 


E ' 

n——oo 




, ( 10 ) 


Equation ([^ for the solute concentration and Eq. Q 
for the flow and pressure fields are solved numerically 
using a boundary integral approach with periodic Green’s 
functions. 

We denote H the fluid domain in a period of the chan¬ 
nel gap, dVl its lower and upper boundaries (the inert and 
active walls), and n the unit normal vector pointing into 
the fluid domain. The two-dimensional periodic Green’s 


where x„ = (a; — ^ + 2 n 7 r, y — t]). These functions may 
be expressed in the closed form 


S^^ = K + ydyK - 1, 

^yy — 7 ^ y^y^^ 

Sxy - y^X^ - Syx^ 

Txxx = “^dxi^K + ydyK), 


T 

xxy 

T 

-Lxyy 

T 

^yyy 

Tijk 


= ‘^dyiydyK), 

= -2ydxyK, 

= 2{dyK-ydyyK), 
^kij '^jkii 


( 11 ) 
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FIG. 2: Dependence of the net flow rate through the channel, 
Q, with the left-right asymmetry of the top wall, 7 , in the 
cases a/H = 1/2, L/H = 27r and 0 < 7 < tt. The four red 
dots correspond to the different panels illustrated in Figurej^ 


for K = i In [2cosh(y) — 2cos(x)]. The computational 
procedure for discretizing the domain boundary is iden¬ 
tical to that used for the diffusion equation 0 . Con¬ 
stant force elements are assumed, singular integrals have 
the singularity removed and computed analytically, and 
non-singular integrals are computed with 16- point Gaus¬ 
sian quadrature.The implementation is based upon the 
authors’ previously published work on the optimal swim¬ 
ming of a sheet [37], with the addition of Tikhonov reg- 
ularisation to improve matrix conditioning. 


III. RESULTS 
A. The role of asymmetry 


(a) Asymmetry y = 0 



(d) Asymmetry y = 3;r/4 



FIG. 3: Solute concentration (colors) and streamlines (lines) 
within the channel with increasing asymmetry in the case 
a/H = 1/2 and L/H = 27r. The four panels above correspond 
to the four red dots in Figure]^ Red (resp. white) streamlines 
correspond to traversing (resp. recirculating) flow regions. 


When inertia and solute advection are negligible, the 
Laplace problem for the solute concentration is linear and 
decouples from the Stokes flow problem, which is also 
linear. Breaking the left-right symmetry is thus required 
in order to create a net flow within the channel, in the 
same way that symmetry breaking is required to achieve 
self-propulsion of autophoretic particles. If the chemical 
properties of the walls are homogeneous, this asymmetry 
can only arise from geometry and therefore, purely sym¬ 
metric profiles such as sinusoidal upper-wall shapes will 
yield zero net flux. 

In order to analyze the effect of asymmetry, we first 
investigate the effect of the asymmetry parameter, 7 , 
on the flow rate. The numerical results are shown in 
Figure At 7 = 0, we recover zero-net flux, as ex¬ 
pected. For asymmetric shapes, we obtain that the flow 
rate within the channel increases monotonically with 7 . 

To gain insight into the origin of this flow rate, we 
illustrate in Figure the dependence of the solute con¬ 
centration distribution and streamlines with 7 . For a 
strictly symmetric profile, 7 = 0 , a flow is induced in 
the channel but one with no net flux. Indeed, a verti¬ 
cal solute concentration gradient is created between the 


two walls due to the fixed-flux emission of solute at the 
upper wall and the constant concentration imposed at 
the lower wall, where the solute is consumed. The upper 
wall is not horizontal, and regions of the upper wall lo¬ 
cated the furthest from the bottom wall are exposed to 
higher concentration than regions corresponding to the 
narrowest channel width. This implies the existence of 
tangential solute gradients along the wall and, hence, of a 
slip velocity that drives a flow within the channel. Due to 
the symmetry of the channel, the flow organizes into two 
counter-rotating flow cells leading to zero fluid transport 
across the channel. For positive activity and mobility of 
the upper wall, the flow is directed away from the active 
wall (i.e. downward) in the regions where the channel 
width is greatest, while it is directed toward the active 
wall (i.e. upward) in the narrowest regions (Figure]^). 

When 7 ^ 0 , asymmetry is introduced in two ways. 
The asymmetric upper wall can now be decomposed into 
a longer backward facing section and a shorter forward¬ 
facing section. The solute gradient on the former is 
weaker than on the latter, leading to a stronger left-to- 
right slip flow along the forward facing section. Addition- 





5 


ally, wall asymmetry increases confinement in the trough 
along the upper wall leading to higher solute concentra¬ 
tion (the rate of production of solute per unit surface is 
fixed). This asymmetry between the wall sections driv¬ 
ing the flow within the channel generates a shape and 
intensity asymmetry between the two recirculation re¬ 
gions, and a traversing streamtube appears (marked by 
dark red streamlines in Figure]^. This streamtube cor¬ 
responds to flow regions that do not recirculate, but are 
transported along the channel, being “pumped” by the 
phoretic mechanism. This tube follows a pattern along 
the channel similar to that of a conveyor belt driven be¬ 
tween the two recirculating regions forced by the slip flow 
on the wall. Along the shorter forward-facing section of 
the upper wall, it is driven by the stronger slip flow that 
dictates the direction of the net flow in this case. The 
tube then separates from the wall where the slip velocity 
changes sign, and circumnavigates around the counter¬ 
rotating flow cell driven by the longer wall section. 

As 7 is increased beyond 7 > 1, the slope of the shorter 
flow-driving wall changes sign, leading to a “folded” ge¬ 
ometry that promotes large confinement effects on the so¬ 
lute concentration distribution (see the difference in color 
scales in Figure]^. This, in turn, enhances the phoretic 
slip and net flow rate. For strong asymmetry, the travers¬ 
ing streamtube is mostly rectilinear and away from the 
active wall, except in a narrowing region where it circles 
around the smaller recirculation region and is driven by 
the phoretic slip within the trough on the boundary. 

This process does not appear to saturate when 7 3> 1 
for fixed amplitude a. In this limit, the flow domain can 
be decomposed into two regions: a complex, thin region 
corresponding to long and thin folds in the wall shape 
where very large concentration gradients are established 
by confinement, and an outer region where a net unidi¬ 
rectional flow is forced within the channel. Beyond ob¬ 
vious practical considerations regarding the manufactur¬ 
ing of such geometries, the assumptions of the current 
model would potentially break down when the asym¬ 
metry parameter 7 becomes too large, as the phoretic 
flow become sufficiently intense for advection to be non- 
negligible (Pe ^ 0). Furthermore, when local concen¬ 
trations become too large, it is likely that the model of 
fixed-flux release would be impacted, and more detailed 
reaction kinetics may be required. 


B. Effect of the pattern amplitude on the flow rate 

The flow within the channel is effectively driven by the 
upper wall, while the no-slip condition on the lower inert 
wall tends to limit the fluid motion. As a consequence, it 
is expected that when the channel gap in the narrowest 
region becomes small (a « H), the net flow rate should 
be small, as the flow viscosity will offer maximum resis¬ 
tance there. However, the corrugation amplitude is an 
essential element to the pumping performance of the de¬ 
vice as it determines the gradient along the upper wall 


(a) Asymmetry y=Tt/A 



Amplitude a/H 

(b) Asymmetry y = ;r/2 



Amplitude a/H 


FIG. 4: The net flow rate through the channel, Q, as a fnnc- 
tion of its relative amplitude a/H in the case L/H = 27r. (a) 
Flux for 7 = 7 r /4 (as Figure!^), showing behaviour represen¬ 
tative of weak asymmetry (unfolded) channels, (b) Flux for 
7 = 7r/2 (as Figure^), showing distinct behaviour for strong 
asymmetry (folded jchannels. 


between the peak and troughs, and therefore the intensity 
of the two recirculating regions driving the flow. When 
a H, it is therefore also expected that the flow rate 
will become negligible. This intuition is confirmed by our 
numerical results in the case of weak asymmetry (7 < 1 , 
unfolded geometry, see Figure]^) for which the net flow 
rate within the channel displays a maximum at inter¬ 
mediate amplitude and decreases to zero in both limits 
a <C Ff and a ~ H. In this case, the limit a H corre¬ 
sponds to a flat upper wall. 

The behavior of the system is however quite different 
when the upper wall is folded (7 > 1 , strong asymmetry, 
see Figure]^). In this case, the limit a H is not limit¬ 
ing to a flat wall, but to a surface with infinitely thin and 
almost horizontal folds. Within these folds, confinement 
creates very large solute concentrations and concentra¬ 
tion gradients. As noted in the previous section, this 
limit is the singular case of a flat wall forced periodically 
by infinitely large slip velocities in infinitely thin regions. 
As a consequence, the net flow rate does not decrease 
to zero for small amplitude, marking a stark difference 
between the folded and unfolded geometries. 
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FIG. 5: Net flow rate through the channel, Q, as a function 
of the channel height H/L for unfolded 7 = 7 r /4 and folded 
7 = tt/2 channels with fixed relative amplitude, a/L = 0.08. 


C. Role of the channel width 

We finally turn to the influence of the third geometric 
characteristic of the channel, namely its mean width-to- 
length ratio. The limit of small minimum width (a « H) 
was already discussed in the previous section and the flow 
rate within the channel vanishes in that limit due to the 
diverging hydrodynamic resistance of the channel. When 
H is large compared to both a and L, the relative concen¬ 
tration distribution along the top wall is not influence by 
the location of the passive wall so the tangential concen¬ 
tration gradients and slip velocity become independent 
of H. As a result the net flow rate through the channel 
varies linearly with the channel height as for a classical 
Couette (shear) flow. This is confirmed by our numerical 
results shown in Figure 


IV. LONG-WAVELENGTH PREDICTION 

When the local height of the channel is small in com¬ 
parison to the typical longitudinal length of the topog¬ 
raphy, i.e. \h'{x)\ 1 , the problem can be solved 

within the framework of lubrication (long-wavelength) 
theory. Defining h{x) = ef{x) with f{x) = 0(1) and 
e= H/L 1, the method consists in solving for the con¬ 
centration and velocity fields as regular series expansions 
in £. On the upper wall, the normal unit vector pointing 
into the fluid is now written 


1 d ‘^1 


£2 dY^ 

1 d'^v 

1 dv 

£ dY dx 

and the boundary conditions at y = ef{x) become 


d'^u 

dx'^ 

d’^v 

dx^ 

du 


dp 
1 dp 

eW’ 


= 0 , 


(13b) 

(13c) 

(13d) 




£ dY 

dx 


+ (14b) 




(14c) 


These, together with the conditions 0 = 0 and u = v = 0 
at Y = 0, suggest searching for solutions of the form 

C{x, Y) = eCi{x, Y) + s^Csix, Y) + 
u{x, Y) = eui{x, Y) + e^u^iyX, Y) + 0{e^), 
v{x, Y) = e^V 2 {x, Y) + 

p{x,Y) = e~^p-i{x,Y) + £pi{x,Y) + 0{e^). 

The flow rate Q is then given by 
rf(^) 

u{x,Y)dY = £'^Q 2 + £‘^Q 4 +0{£^), (16) 


(15a) 

(15b) 

(15c) 

(15d) 


Q = £ 

with 


Qj = [ Uj-i{x,Y)dY. 
Jo 


(17) 


The flow is incompressible and steady, therefore Q and 

Qj do not depen d on x. _ _ 

Inserting Eq. (15a) into Eqs. (13a I and (14a I gives at 
leading order 

Ci(x,Y) = Y. (18) 

Eqs. ( 13b I, ( 14b ) and ( |15b[ ) then provide at 0 (£): 

U4{x,Y) = ?^{Y^-Yf) + ^, (19) 

with p_i the leading-order pressure distribution which 
is vertically invariant. The function p-i{x) is periodic, 
therefore 


Q 2 = 0 and p-i{x) =-6/f. 


( 20 ) 


~ ~7TTW~ 

p2 T 2 p.3 f 3 

= —By + £f'ex H-^—Gy- ^—Bx Y O(e^). (12) 


We see that in the lubrication limit, a velocity field is 
present at 0 (£), which takes the form of two recirculating 
regions, but does not give rise to any net flow through the 
channel at this order. After substitution and application 
of the continuity equation, we obtain 


Defining a rescaled vertical coordinate y = £y, the 
Laplace and Stokes flow problems are now given by 


1 d'^C d'^C 


£2 dx^ 


= 0 , 


(13a) 


wi(a;,y) = 3(^|^)r2-2(^^^y, ( 21 a) 


( 21 b) 
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tion (long-wavelength) theory predictions for 7 = tt/IO and 
with pi(x) a function of x only. Integrating the previous a = H/3 as & function of an increasing slenderness, 1/e. 
equation in Y finally provides the flow rate at 0{e^) 


Qi = 2 


3/5 

40 







12 " 


(24) 


Using the periodicity of pi, Q 4 can be computed by di¬ 
viding the previous equation by p, taking the spatial 
average in x, and integration by parts. The result can be 
rewritten in terms of the original channel height h{x). At 
leading order we obtain that the flow through the channel 
is given by 


Q = 


11 




30 (l//i3) 


(25) 


where (•) is the spatial average over a period. 

We see two important results: (i) the flow rate is in¬ 
timately linked to the distribution of local slope along 
the wall h'{x), and (ii) shape asymmetry is essential for 
the existence of a net flow. Indeed, slip flow along the 
active wall arises from the orientation of the wall with a 
component along the leading order solute concentration 
gradient. A non zero h'(x) is therefore sufficient to guar¬ 
antee the existence of a local flow but not necessarily of 
a net flow through the channel. This separation of scales 
is clearly visible in the lubrication expansion: the leading 
order flow arises from the local channel geometry (i.e. the 
fact that the wall is not flat and orthogonal to the lead¬ 
ing order solute gradient). However, at this order, the net 
flow is zero because the flows driven by the forward- and 
backward-facing walls exactly cancel out. A net flux re¬ 
sults from an imbalance between these local flows which 
can only be induced by geometric asymmetry. Left-right 
symmetric profiles are characterized by an even channel 
height, h{x) = h{—x). Consequently the function h /h? 
is odd and thus exactly averages to zero, so that (^4 is 
identically zero (all higher orders are expected to be zero 
as well). 


The result in Eq. (25) is a weighted algebraic spatial 
average of the slope of the active wall. More precisely. 


the leading order flow rate through the channel, Eq. (251, 
is the ratio of two integrals; the numerator is the mean 
flow forcing due to the asymmetry of the channel, while 
the denominator is the average hydrodynamic resistance 
of the channel over a wavelength. 

This leading-order prediction is compared to the full 
numerical simulations in Eigure Eor a fixed asym¬ 
metry parameter 7 and relative amplitude a/iJ, several 
simulations are performed for increasing LjH (note that 
as H is reduced, a is reduced in the same amount) and 
the flow rate through the channel is shown to converge 
for large slenderness to the prediction of the lubrication 
theory. 

Note that the lubrication result, Eq. (25), is valid for 
any ratio iL/L I, regardless of the relative magnitude 
of the mean channel height H and the perturbation am¬ 
plitude a. In the limit of small wall roughness (a ^ H), 
the hydrodynamic resistance (the denominator in Eq. 25) 
is independent of a at leading order and simply scales 
as 1/H ^, while the phoretic forcing (the numerator in 
Eq. scales as As a consequence, Q scales as 

a^H when a ^ H. 

In the opposite limit a H, using classical asymp¬ 
totic expansions to compute the leading order contribu¬ 
tion to the integrals in Eq. (25) [38] , one can show that 
the flow forcing due to the channel’s asymmetry (i.e. the 
numerator in Eq. 25) remains finite and 0{H), but that 
in contrast the hydrodynamic resistance diverges. More 
specifically, a standard lubrication calculation leads at 
leading order to 


I \ 3v^ 

P/ “ I6iLi/2(i7_a)5/2’ 


(26) 


As a consequence, Q ~ — a)^/^ when {H — a) 

H. 

Note that since the limiting factor is then the channel 
width at the narrowest point, and its impact on the hy¬ 
drodynamic resistance, it is expected that this scaling in 
{H — a)5/2 should hold true even when H is not small 

























and could be recovered through a new lubrication expan¬ 
sion limited to the narrow-gap region of the channel |39j , 
provided the curvature of the wall in that region remains 
finite. 


clockwise flow through a two-dimensional annular chan¬ 
nel driven by the geometric asymmetry of the inner active 
wall. 


V. CONCLUSIONS 
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particles has demonstrated that fuel-based mechanisms 
represent a promising route to designing self-propelled 
systems that rely only on chemical reactions and the in¬ 
teraction with the immediate environment to create lo¬ 
comotion. The results presented in this work show that 
this is also true for the dual problem of pumping flow 
within a micro-channel, and that geometric asymmetry, 
rather than chemical patterning of the channel walls, is 
sufficient to create a net flow. Our results provide insight 
into the flow dynamics within the channel, and the mech¬ 
anism leading to the net fluid transport: the breaking of 
symmetry between two recirculating flow regions driven 
by wall slip velocity, and the emergence of a conveyor- 
belt-like flow within the channel. 

For simplicity, we focused in this paper on a reduced 
set of wall shapes with one active wall and the other 
one passive. The numerical methodology and the long- 
wavelength theory, Eq. (25), are however valid for any pe¬ 
riodic channel geometry in two dimensions. Furthermore, 
these results were obtained within the simplified frame¬ 
work of a fixed-rate release of solute by the active wall. 
Previous studies on autophoretic self-propelled particles 
have shown that the exact reaction kinetics, in particular 
the dependence of the reaction rate on the local solute 
concentration, may significantly impact the system dy¬ 
namics [S3I1S]. We expect for example the direction of 
pumping to be impacted by reaction kinetics, although 
the basic result showing the emergence of a net flow due 
to geometric asymmetry of the phoretic wall should re¬ 
main true. Finally, our study focused on the particular 
case of self-diffusiophoresis. Because of the formal simi¬ 
larities in the equations of the problem, these results can 
be generalized easily to other phoretic mechanisms such 
as self-thermophoresis or self-electrophoresis. 

Looking forward, the results of this work could be gen¬ 
eralized to a larger range of geometries, including closed- 
loop channels for which pressure-driven flows can not 
easily be achieved. This is shown in Figure where we 
have adapted our numerical approach to compute the net 


SM acknowledges the support of the French Ministry 
of Defense (DGA). TDMJ is supported by a Royal Com¬ 
mission for the Exhibition of 1851 Research Fellowship. 



0 0.5 1 1.5 2 2.5 3 3.5 

FIG. 7: Concentration distribution and flow streamlines 
within an annular closed-loop channel with a geometrically- 
asymmetric inner active wall releasing solute with a fixed flux, 
and a passive circular outer wall with uniform concentration. 
The recirculating streamlines are shown in white while the 
traversing streamlines are plotted in red, and correspond to a 
clockwise-rotating flow. 
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